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Abstract 

We present a model for thin film growth by particle deposition that takes 
into account the possible evaporation of the particles deposited on the surface. 
Our model focuses on the formation of two-dimensional structures. We find 
that the presence of evaporation can dramatically affect the growth kinetics 
of the film, and can give rise to regimes characterized by different "growth" 
exponents and island size distributions. Our results are obtained by extensive 
computer simulations as well as through a simple scaling approach and the 
analysis of rate equations describing the system. We carefully discuss the 
relationship of our model with previous studies of the same physical situation, 
and we show that our analysis is more general. 
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Typeset using REVT^X 



The growth of nanostructures and thin films prepared by atomic deposition is recognized 
to be both a standard and a promising way to prepare new materials |Il-f||- Several models 
|5|-|| have led to a good understanding of the growth properties in the simplest cases, when 
only a limited number of physical ingredients are included in the simulations: Deposition, 
Diffusion and Aggregation as in the DDA model ||. A drawback of this conceptual sim- 
plicity is that the range of experimental situations accurately described by these models 
is limited. It is actually limited to (beautifully) artificial experimental setups where great 
care is taken to avoid complications (contamination, surface defects, etc. see for example 
|T0|-|T2|). Clearly, many (technologically) interesting experimental situations are much more 
complex. Progress towards their understanding demands the inclusion of processes that have 
been left out of the first models. For example, including reversible aggregation |13| allows 



to understand the saturation of island density before island-island coalescence and produces 
compact islands. In this paper we show the effects of including evaporation of the atoms 
from the surface. 

Evaporation, i.e. the possibility of desorption of adatoms from the surface, is a feature 
that should be observed for any system at high enough temperatures. In this sense, it is a 
phenomenon that is as general as the rest of the ingredients of recent models of film growth, 
and is capable of completely changing the quantitative behaviour of the system (scaling 
dependence of island density on the deposition parameters, island size distributions ...). 
Moreover, the effects of evaporation have already been studied experimentally [pL4|— p76 1 . Thin 



film growth models which include evaporation have already been studied using a mathemat- 
ical analysis of rate equations P, p!5lJl7| -[19| . Computer simulations of such models, aiming a 



quantitative analysis of scaling relations or island size distributions have, to our knowledge, 
never been carried out. The point is that computer simulations represent an "exact" way 
of reproducing the growth, in the sense that they avoid the mean-field approximations of 
mathematical rate-equations approaches P J20|J21[] . We find that a careful consideration of 
the downward transport of monomers from islands leads to results that significantly differ 
from previous studies [|l5| ,|r^-[l9l. 



Evaporation is of course not the only process that needs to be added in order to better 
describe most of the film growth experiments. Other candidates include interlayer transport, 
particle dissociation (if the incident particles are molecules or clusters rather than atoms) 
or intricate chemical interactions between adatoms and the surface (for a taste of the real 
growth world, see for example refs. J2"2"| , |2"3" |). These effects are less pervasive than evaporation, 
but nevertheless, should eventually be considered, keyed to specific experimental systems. 

The paper is organized as follows. Section I briefly presents the model and discusses 
some of its approximations. Then, in section II, we give a rapid overview of the growth of 
the films when evaporation is taken into account, and try to give a physical intuition for 
what is going on. After this, section III presents a simple scaling approach which gives the 
behaviour of the saturation island density as a function of the parameters. This scaling 
approach is completed in section IV by a more rigorous approach based on a careful analysis 
of the rate equations of the system. Section V confirms the preceding results by extensive 
computer simulations of the model. Finally, in section VI, we discuss our results and the 
differences with the previous approaches [P75],P~T^P~9[ . 



I. PRESENTATION OF THE MODEL 

In this work we will describe the properties of a still oversimplified submonolayer thin 
film growth model which includes the four most important physical ingredients of these 
systems: 

(1) Deposition. We will assume that atoms are deposited at randomly-chosen positions 
of the surface at a flux F per unit surface per unit time. 

(2) Diffusion. Isolated adatoms can move in a random direction by one diameter, or one 
lattice spacing, which we will take as our unit length. We denote by r the characteristic 
time between diffusion steps. 

(3) Evaporation. Isolated adatoms can evaporate off the surface at a constant rate. We 
denote by r e the mean lifetime of a free adatom on the surface. As an approximation, we 
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will assume that the these desorbed atoms do not return to the surface. It is useful to define 
X s = \[Dr e the mean diffusion length on the substrate before desorption. 

(4) Aggregation. If two adatoms come to occupy neighboring sites, they stick to form 
an island. We will mainly consider the case of irreversible aggregation, but a more general 
analysis of reversible aggregation, assuming for simplicity the existence of a critical size i* 
below which islands tend to dissociate is given. In any case, islands are assumed to be 
immobile and do not evaporate. 

In the following, we call particles or adatoms the isolated atoms (or monomers) that 
are deposited on the surface, and islands a set of connected particles (thus excluding the 
monomers) . 

Some remarks on the assumptions of this simple model regarding its connection to the 
experiments are now addressed. 

(a) Second layer — When a particle "falls" on top of an island, we assume that the particle 
deposited on the second layer moves and evaporates essentially as any other particle, except 
that its diffusion constant and evaporation rate now correspond to the process ocurring on 
a substrate of the same element as the deposited particles. Practically, the mean diffusion 



length on the island before desorption, defined as X s = y D*r*, where the * indicate the 
values on the island, can be different from Xs- As the simplest scenario, we assume that 
when the particle reaches the border of the island, it immediately jumps down and increases 
the area of the island. We discuss when this effect can be ignored without affecting the 
scaling results, which leads to applications to situations in which downward transport of 
monomers from islands is highly improbable. This could be the case when: (i) there is 



a barrier [24[] at the edges of the (first layer) clusters which prevents single particles from 
"falling" on the substrate and/or (ii) Particle diffusion on the second layer is much smaller 
than diffusion on the substrate. 

(b) Island diffusion — We neglect in this model the possibility for dimers, trimers or larger 
islands to diffuse on the substrate. Island diffusion has been observed in experiments [^7j] and 



molecular-dynamics simulations |28|]. The effects of island diffusion on the growth of thin 



films without evaporation have been addressed in Refs. [fl, p0| , pl] ,j5 p9| , |30| . While these effects 
are significant since they change the growth exponents, they do not modify fundamentally 
the growth mechanisms. We assume that the same is true when evaporation is included. 

We study here the first stages of the growth, roughly until the number of islands on 
the substrate saturates. The reason is that it is in this stage that such a simplified model 
can be of some help to experimentalists who want to understand the microscopic processes 
present in their experiments. These fundamental microscopic processes are most easily 
detected in the first stages of the growth, since in the subsequent stages additional processes 
can be involved (additional diffusion paths, interlayer transport, geometrical details of the 
lattice. . . ). 



II. QUALITATIVE DESCRIPTION 

Before going into the details of the calculations and their confirmation by computer 
simulations, we present a simple picture of the growth mechanisms of the submonolayer 
structures under consideration. For simplicity, we only consider here the case i* = 1, i.e. 
irreversible aggregation and we neglect the particles deposited on top of the islands (despite 
the apparent lack of generality of this last hypothesis, the results obtained here are quite 
general, as we show in section IV). 

The qualitative evolution of the system is essentially as follows. The system initiates 
as a clean empty surface. Monomers are then deposited at a constant rate on the surface 
and are allowed to diffuse and/or desorb (evaporate). When two (or more) monomers meet, 
they aggregate irreversibly to form a static island. As more of these encounters occur, the 
island density increases with time and islands become larger by capturing adatoms. At some 
point, islands are so large that they start touching (coalescing), and monomers are rapidly 
captured. These two effects lead to a saturation in the number of islands. Interestingly, the 
saturation is attained when the surface coverage reaches a value close to .15, independently 
of the parameter values. This last point is discussed in detail below. We now turn to a 
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more detailed discussion of the evolution of the systems in two limiting cases : complete 
condensation (evaporation is negligible) and strong evaporation. 

First consider the situation where evaporation is negligible. This means that atoms 
deposited on the surface almost never evaporate before aggregating (after this, they are safe 
since islands do not evaporate). This situation can be expected to happen when X s ^> £ 



where Xs = y/Dr e is the adatom diffusion length on the substrate before desorption and 
i is the typical distance between islands. Fig |l|a shows the evolution of the monomer 
and island densities as a function of deposition time. We see that the monomer density 
rapidly grows, leading to a rapid increase of island density by monomer- monomer encounter 
on the surface. This goes on until the islands occupy a significant fraction of the surface, 
roughly 1%. Then, islands capture rapidly the monomers, whose density decreases. As a 
consequence, it becomes less probable to create more islands, and we see that their number 
increases more slowly. When the coverage reaches a value close to 15%, coalescence will 
start to decrease the number of islands. The maximum number of islands N max is thus 
reached for coverages around 15%. Concerning the dependence of N max as a function of the 
model parameters, the DDA and related models have shown that the maximum number of 
islands per unit area formed on the surface scales as N max ~ (F/D) 1 ^ 3 , or £ C c — (F/D)^ 1 ^ 6 
|5| Jl8|Jl9t| where CC stands for Complete Condensation. These values are independent of r e 
since evaporation is not significant. 
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FIGURES 




FIG. 1. Evolution of the monomer and island densities as a function of the deposition time (in 
monolayers) for different limits : (a) complete condensation, F = 10~ 8 , r e = 10 10 (r = 1). These 
values mean : Xs = 10 5 and ice = 22 (b) important evaporation, F = 10 -8 , r e = 600 (r = 1). 
These values mean : Xs = 25 and ice = 22. In (b) the "sticking" curve represents the total number 
of particles actually present on the surface (the coverage) divided by the total number of particles 
sent on the surface (Ft). It would be 1 for the complete condensation case, neglecting the monomers 
that are deposited on top of the islands. The island density curve has been shifted vertically by 
(+1). The solid line represents the constant value expected for the monomer concentration, while 
the dashed line corresponds to the linear increase of the island density (see text). 

When r e decreases (i.e. as the evaporation rate increases), Xs decreases and eventually 
becomes smaller than ice- in this regime, evaporation is going to significantly alter the 
growth dynamics, as shown in Fig [L|b. The main point is that now the monomer density 
becomes roughly a constant, since it is now mainly determined by the balancing of deposition 
and evaporation. As expected, the constant concentration equals Fr e , as shown by the solid 
line. Then the number of islands increases linearly with time (the island creation rate is 
roughly proportional to the square monomer concentration). We also notice that only a 
small fraction (1/100) of the monomers do effectively remain on the substrate, as shown 
by the low sticking coefficient value at early times (the sticking coefficient is the ratio of 
particles on the substrate (the coverage) over the the total number of particles sent on the 
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surface (Ft)). This can be understood by noting that the islands grow by capturing only the 
monomers that are deposited within their "capture zone" (comprised between two circles of 
radius R and R + Xs if we neglect X s , see Figure 0). The other monomers evaporate before 
reaching the islands. As in the case of complete condensation, when the islands occupy a 
significant fraction of the surface, they capture rapidly the monomers. This has two effects 
: the monomer density starts to decrease, and the sticking coefficient starts to increase. 
Shortly after, the island density saturates and starts to decrease because of island-island 
coalescence. 




FIG. 2. Schematic capture zone for an island in presence of atom evaporation. X s stands 
for the mean length before desorption for an atom diffusing on top of an island, whereas Xs 
corresponds to the same parameter on the substrate (see text). The sizes of the capture zones will 
be justified below, in section IV. 

If r e is further decreased (precisely r e < r i.e. the time a particle remains on the surface 
is less than the time it needs to move), then, clearly, diffusion plays no role ("Xs < 1")- 111 
this situation, islands are formed by direct impingement of incident atoms as first neighbors 
of adatoms, and grow by direct impingement of adatoms on the island boundary. This 
situation, although apparently uncommon, is not physically impossible and it also allows us 
to test our predictions over a larger range of parameters. 
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III. SCALING ARGUMENTS 



In this section we present simple scaling arguments that allow to find the dependence 
of the maximum island density N max as a function of the deposition parameters (Flux F, 
Diffusion time r and Evaporation time r e ). These arguments were originally formulated 
in || for the special case of deposition without evaporation on a high-symmetry terrace. 
Here, the argument is extended to the case of non-negligible evaporation. We recall that the 
atomic size is taken as the length unit. For simplicity, we neglect in this section the effects 
of deposition on the second layer : they will be studied in great detail in the next section. 
We will show there that the regime studied in this section corresponds to a large range of 
physical situations. 

The first stage of the argument requires the determination of the nucleation rate per 
unit surface and time, l/r nuc . A nucleation event takes place when two adatoms meet. This 
happens with a probability per unit time Dp 2 , where D = l/(4r) is the adatom diffusion 
constant, and p the adatom density. Thus, 

— w Dp 2 (1) 

Time 

Another, independent equation can be written down to relate the nucleation rate and 
the island density N. It states that in the area I 2 = 1/N occupied by an island, only one 
(on average) nucleation event takes place, during the time t c needed for the growing islands 
to come into contact. Thus, 

1 N (2) 



Time t c 

The time t c is readily computed by knowing the growth velocity of an island, which 
in turn requires the knowledge of the adatom density. We consider in the following the 
situations of interest for this paper. 
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A. No evaporation, compact 2-D islands 



The adatom density results in this case from a balance between deposition at a rate F 
and capture by the stable islands at a rate DpN, so that 

p w F/(DN) pa F^/D (3) 

The area of an island of linear size R increases by capture of all adatoms falling in an 
area £ 2 , so that R 2 pa F£ 2 t. At t = t c , R w 4, and thus £ c » Using (|) and (g) one 

finds Dp 2 « FN. From (0), JV « F£*/D, or 



' (4) 

The value 1/3 corresponds to the special case i* — 1 of the general formula i*/(2 + z*) 
for the exponent ||[]]|[l9|. 



B. Strong evaporation 

Strong evaporation means the adatoms are more likely to disappear due to desorption, 
with probability l/r e per unit time and site, than to be captured by an island. Therefore, 
the adatom diffusion length before desorption, X$ = ^Dr e , is shorter than the average 
island-island distance, £ s . In this case, the adatom density results from a balance between 
deposition and desorption at a rate p/r e , so that 

p^Fr e ^ FX 2 S /D. (5) 

A growing island is only able to capture the adatoms falling at a distance smaller than 
the adatom diffusion length before desorption, X$ (see Figure 0). If X$ is smaller than the 
island size R, the island will capture the adatoms falling inside an annulus of width pa X$ 
around its border, as well as those directly impinging on its edge. Since the area of the 
annulus is pa RX S if R > X s , at time t one has R pa F(l + X s )t. At t = t c , R pa £ s , and 
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thus t c ~ £ S /[F(1 + Xs)] ~ l/[-FA^ 1 / 2 (l + Xs)]- Direct impingement is also important for 
nucleation, so that (D) becomes 

(F + Dp)p « F(l + X 2 s )p w F 2 r e (l + X s ) 2 (6) 

where we also used @. 

From (0) and (§), one finds F 2 r e (l + X s ) 2 w F(l + X 5 )iV 3 / 2 , or 

iV « (Fr e ) 2/3 (1 + X 5 ) 2/3 . (7) 

For a comparison with previous results, it may be worth it to rederive (0) for any value 
of %* . The nucleation rate reads, for any i*, {F + Dp)N i * ) where iVj* is the density of critical 
nuclei (clusters of size i*). Following Refs. jl8| , [l9fl , we assume that Np satisfies Walton's 
relation Ni* ~ p e . Inserting into @ yields l/r nuc w F(Fr e )**(l +X S ) 2 . Finally, using @ 
yields F(Fr e f(l + X s f w F(l + X 5 )iV 3 / 2 , or 

iV^ (Fr e ) 2l * /3 (1+X s ) 2/3 . (8) 

Note that the validity of Walton's relation in the presence of strong desorption is highly 
hypothetical. It is a point which we reserve for future work. In the rest of the paper, we 
only consider the case i* = 1. In section IV, the results for i* — 1 are derived by a more 
rigorous rate-equations approach. 

C. Crossover scaling 

We will derive here an approximate analytic expression for the crossover scaling func- 
tion connecting the complete condensation, diffusion and direct impingement regimes de- 
scribed previously. In other words, we will compute the maximum island density N max as 
a function of F, r and r e . We will then show that, if we measure all lengths in terms 
of Xs, N max (F,T,T e ) is a function of its arguments only through a special combination: 
XsN max = g e va P (XsN evap ) } where N evap = [(1 + X 5 )Fr e ] 2/3 , and g e va P (x) satisfies 
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X 1 / 2 f or x _> 00) 

9evap(x) rsj <j (9) 

x for x — > 0. 

To do this, we assume that the islands are large enough with respect to the atomic size, 
that we can neglect the curvature of their boundary. Then, the adatom density in the region 
between two islands whose edges are at a distance £, obeys the equation 



p = F + DV 2 p - - 

r. 



(10) 



with the boundary conditions p(±£/2) = (the origin being midway between the islands). 

In the quasi-stationary approximation, p m 0, and equation ( [TDD can be solved. The 
solution reads 



p(x) 



Ft 



1 - 



cosh(ra) 



(11) 



cosh(K£/2)_ 
where k = 1/Xg 

This formula is needed to compute the nucleation rate Dp 2 . The latter is a mean field 
quantity, independent of x. Letting thus x = in (|TT|) -since the highest nucleation proba- 
bility is at the terrace centre, given the symmetry of our problem-one finds 



p = Ft 



1 



cos1i(k£/2) 



2Ft 



sinh 2 (K£/4) 
cosh(/t£/2) 



(12) 



where we used the identity cosh(x) — 1 = 2 sinh 2 (x/2). 

The next task is the determination of the total island density N. Its time variation is 
simple: N increases each time a new island is nucleated, so that (N)i = Dp 2 . On the other 
hand, N decreases when two islands touch and coalesce. Following previous authors ||TT|,|T§ 



we write (N) 2 = —(dA/dt)N 2 , where A ~ R 2 is the average area of an island of linear size 
R. This means that coalescence results from binary encounters of immobile islands, whose 
area increases at a rate dA/dt. Collecting (N)i and (iV) 2 yields 



N = Dp 2 — ^N 2 
p dt 



(13) 
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By definition, N max , which is what we are interested in, satisfies N = 0. This means that 
the maximum of the island density can be found by balancing the nucleation rate against 
the coalescence rate, and that the balance is reached when the island size R ~ £. One can 
thus write 

Dpi "&) B J l - (14) 

The growth rate of an island can be computed assuming that its radius is large with 
respect to the atomic size. One can then treat the island edges as straight steps, which 
yields 

— « FX s tanh(/^/2) . (15) 

(JjL 



Since dA/dt ~ RdR/dt, one finds 

dA, 



\R= 



FX s £tanh(Ki/2) . (16) 



dt 

Note that the two limits (dA/dt)\ R=l » Fl 2 and (dA/dt)\ R=e « F£X S are well repro- 
duced for < 1 and for k£ 3> 1, respectively. In fact, in the limit X$ — > 0, the area 
increase is FR, due to direct impingement. In order to interpolate all the way through to 
this regime, we finally write 

dA 



\R= 



« F{1 + X s )£tanh(Ki/2) . (17) 



dt 

In the limit where the diffusion length X$ becomes very small, the nucleation probability 
is also determined by direct impingement of a beam atom on a nucleus, which is described 
by F p. Then, the nucleation term in (TA) reads 



(F + Dp)p w F 



l+X s 



2 sinh 2 (K£/4) 



Ft, 



smh 2 (K£/A) 
cos1i(k£/2) 



(18) 



cosh(fc£/2) 

In fact, a simpler approximate expression will be used, which reproduces the correct limiting 
behaviour at large and small X$: 



[l + X s f L sinh 2 (/^/4) 



(F + Dp)p 



F T 

e cosh(K£/2) 



(19) 
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Finally, (0) and (0) yield 



Ft, 



sinh 2 (/^/4) 
cosh(/t£/2) 



F(l + X s )£ta,uh(Ki/2)N, 



(20) 



or, letting i = N'^ 2 , 

(l + X s )(Fr e 



cos 



h(KiV-V 2 /2) 



tanh(«:Ar-V 2 /2)iV3/ a 2 r . 



(21) 



Lsinh^iV-^V^J 

This is the announced crossover scaling formula, in implicit form. It can be cast in the 
form: 



[l + X s )Xl(Fr e 



cosh(/dV-i/ 2 /2) 



-i 2 



sinh 2 (KAT-i/ 2 /4)J 



tanh( K iV-V 2 /2)(X|iV max ) 3 / 2 



or, 



;i + x 5 )x|(^-j =f euap (x s N£ e 



where 



fevap \%) 



cosh (i) 



tanh | — | 



Letting N evap = [(1 + X s )FT e } 2 ^ 3 , and inverting /, one finds as promised: 



(22) 



(23) 



(24) 



XsN max — g eva p{XgN ( 



evap ) 



(25) 



IV. RATE EQUATIONS 

In this section we will study the model taking into account the particles deposited on 
top of the islands with their specific diffusion and evaporation parameters. We will show 
that for a large range of parameters, the exponents are the same as those predicted by the 
preceding analysis. We use a "rate-equations" approach which has been shown to give a 
good description of the submonolayer regime |G|,|£^ |r7|,|I§||5I| . For the sake of completeness, 
we include a detailed calculation of the "cross sections" for monomer capture. These "cross 
sections" have also been calculated in [p,15|,17| for example. 
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We will not keep careful track of many of the numerical geometric constants. The rate 
equation describing the time evolution of the density p of monomers on the surface will be, 
to lowest relevant orders in F: 

^ = F (i _ o) _ L - F p - 2a oP - a t N (26) 
at r e 

The first term on the right hand size denotes the flux of monomers onto the island free 
surface, (9 is the island coverage discussed below). The second term represents the effect 
of evaporation, i.e. monomers evaporate after an average time r e . The third term is due to 
the possibility of losing monomers by effect of direct impingement of a deposited monomers 
right beside a monomer still on the surface to form an island. As discussed above, this 
"direct impingement" term is usually negligible, and indeed will turn out to be very small 
in this particular equation, but the effect of direct impingement plays a crucial role in the 
kinetics of the system in the high evaporation regimes. The last two terms represent the 
loss of monomers by aggregation with other monomers and with islands respectively. The 
factors a and <7j are the "cross sections" for encounters and will be calculated below. 
The number N of islands will be given by: 

f = Fp + .„p (27) 

where the first term represents the formation of islands due to direct impingement of de- 
posited monomers next to monomers already on the surface, and the second term accounts 
for the formation of islands by the encounter of monomers diffusing on the surface. 

For the island coverage 9 i.e. the area covered by all the islands per unit area, we have: 

d9 

— = 2[Fp + a oP ] + a t N + JN (28) 

The term in brackets represents the increase of coverage due to formation of islands of 
size 2 (i.e. formed by two monomers) either by direct impingement or by monomer-monomer 
aggregation. The next term gives the increase of coverage due to the growth of the islands 
as a result of monomers aggregating onto them by diffusion, and the last term represents the 
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growth of the islands due to direct impingement of deposited monomers onto their boundary, 
or directly on the island. The estimation of the value of J is carried out below. The total 
surface coverage is given by 9 + p ~ 9 except at very short times. 

The next step in the analysis consists in estimating the cross sections a Q and crj. This 
is done by evaluating the diffusive flux of monomers into a single single solitary island 
represented by an absorbing disk of radius R centered at the origin (throughout this analysis 
R stands for the typical island radius). The typical radius R of the islands will be taken to 
be: 

1/2 

N 



The cross sections are evaluated in the quasistatic approximation, which consists in 
assuming that R does not vary in time and that the system is at a steady state. Thus, we 
have to solve the equation 

fd 2 P ldP\ 1 

( dr 2 ^ r dr J ^ r e ^ ^ ^ 

subject to the boundary condition P(r = R) = 0, where -P(r) is the monomer concentration 
at position r and D is the diffusion constant of the monomers. The solution is given by: 

where X$ = (-Dr e ) 1//2 can be identified as the characteristic distance traveled by a monomer 
before evaporation, and Kq(x) is the modified Bessel function of order zero. We can rewrite 
Eq. |3l] as: 

P(T) = " I 1 " IwvW)) ' (32) 

where p is the monomer density appearing in the rate equations, so that P(r) — > p far from 
the islands. This, again, is a mean field way of including the effects of the other islands 
in the flux of monomers into the island in consideration. The cross section Oi can then be 
calculated as the total diffusive flux of monomers into the boundary, 
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Oi = ^BD\ — \ =2*Dp( — ) * ' (33) 



drJ r=R ~ ' \X S ) K (R/X S ) 

The cross section for monomer-monomer encounters a Q is obtained from the same formula 
substituting R by the monomer radius, and D by 2D as corresponds to relative diffusion. 

The rate equations, with the explicit expressions for <7j and a are still too complicated 
to solve exactly. At best we can focus on extreme cases. The first step then is to evaluate 
the extreme values of U{ for the case in which R Xs and R <C Xs- Using the known 
asymptotic values of the Bessel functions, and once again omitting numerical constants, we 
have 

' %2 for R » X s , 
Oi ~ | (34) 

>m§Sm {orR«X s . 

From here on we will neglect the logarithmic variations appearing in the cross section. To 
evaluate a we notice that the monomer radius is a small constant (it is usually taken as 
half the unit of length), so for all values of Xs, to a good approximation we have 

°o ~ Dp (35) 

Formally, the diffusive cross section for monomers should vary in the case Xs <C 1, but 
this is of no consequence since in that regime the contributions from diffusion are negligible. 

Now we turn to the evaluation of the direct impingement flux J. As was the case above, 
the exact evaluation of J is rather difficult, as it results from the solution of a moving 
boundary problem. Thus, once again, we must resort to a quasistatic approximation. 

The flux J will consist of two terms, Jb the contribution due to direct impingement 
at the exterior boundary of the island, i.e. on the substrate, and Jj the contribution due 
to particles falling on the island, diffusing to the boundary, and "falling" over the edge to 
increase the size of the island. We will continue to assume that the islands are compact, thus, 
the contribution due to direct impingement at the boundary will be given by FR (where we 
have omitted geometrical constants). To evaluate the second contribution we must solve: 
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D* 



d 2 P 1 dP 



) 



+ F- — P(r) = 



(36) 



subject to the boundary condition P{r = R) = 0. We recall that L>* and t* represent the 
values of the diffusion constant and evaporation time of the particles on a substrate of their 
same species (clearly, for homoepitaxy D* = D and r* = r e ). The boundary condition 
corresponds to the simplest situation, in which there is no barrier at the edge of the island 
and every particle reaching the boundary falls. The increase in coverage will be given by the 
total flux accross the barrier, i.e. Jj ~ —RD*^. In writing this equation we have assumed 
that no nucleation occurs on top of the island, for this to be the case, we must require that 
Pmax < 1 on the island, which will be the case for small enough F. 

The solution to the above equation can be readily found, and from it we find that 



where X* = ^ D*r* is the typical distance a monomer can diffuse on an island before 
desorption, and Iq(x) and I\{x) are modified Bessel functions. Again using the known 
properties of Bessel functions, we can distinguish two limiting behaviors: 



We are now in a position to analyse the limiting cases of the system described by the 
rate equations. First we consider the new behaviour brought about by the presence of 
evaporation. We define the high evaporation regimes as the systems in which all but the 
first two terms on the right-hand side of Eq. ^6| are negligible. Then, since the coverage is 
small, we have 




(37) 





(38) 



A. High Evaporation Regimes 




(39) 
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and p reaches a steady state value p = Fr e in a time of order r e . Thus, as mentioned 
above, the high evaporation regimes are characterized by having a constant monomer density 
(after an initial transient) throughout most of the evolution of the system. Under these 
circumstances Eq. |27] becomes trivial and predicts that the number of islands at time t is 
given by 



N ~ F 2 r P 



1 + X 



t (40) 



If the evaporation rate is very high, so that X$ <C 1, then even the smallest islands at 
the earliest stages of evolution will satisfy the relation R ^> Xs (will refer to this situation 
as extreme incomplete condensation). If we assume that the diffusion length on the islands 
X* is also much less than one, then Eq. |28| becomes 

~ FRN ~ F{9Nfl\ (41) 

etc 

where we have neglected the increase of coverage due to the formation of small new islands 
in Eq. |28|, as most of the coverage is due to the large islands. And the number of islands 
evolves as 



N ~ F 2 T e t, (42) 



Solving for 6, we obtain 



9 ~ F 4 r e t 3 . (43) 



Now comes a crucial assumption, namely that coalescence occurs when the coverage reaches 
a fixed constant value, say 6 C = 1/6, and that this value is essentially independent of the 
details of how the surface is covered. Thus, the time t c at which coalescence occurs will be 



t 



(FVe)" 1/3 (44) 



Since the number of islands is a monotonically increasing function of time up to the time of 
coalescence, we can estimate the maximum number of islands on the surface as 

19 



N max ~ N(t c ) ~ (Fr e ) 2/3 (45) 

If instead we consider the case in which i? 3> X* 3> 1, and X s <C 1, then island formation 
will still only occur via direct impingement, but the increase of coverage will be primarily 
due to monomers landing in the capture zone on the islands. In this situation, we will have 

^ ~ X*FRN, (46) 



and again 



N ~ F 2 r e t, (47) 



from which we obtain 



6 ~ F\ e X*t\ (48) 

The coalescence time can be estimated as before, leading to a maximum number of islands 
that scales as 

N max ~ iV(t c ) ~ (Fr e X;) 2/3 (49) 

If we instead consider the case of extreme mismatch between the materials, so Xs <C 1 
but X s ^> R throughout the evolution of the system, the the capture zone on the islands 
becomes the whole island, at which point the coverage increases as 

f~F», (50) 

where we used NR 2 ~ 9. Coalescence then occurs at at a time t c ~ 1/F, and the maximum 
number of islands at that time will be 

Nmax ~Fr e , (51) 

Different and more interesting regimes exist for the situation in which Xs ^> 1, still in 
the high evaporation regime, i.e. the evolution of the monomer density is still essentially 
controlled solely by the deposition and evaporation. This requirement actually limits the 
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value of Xs- The upper bound for the range of Xs in which this situation prevails is discussed 
below. First we consider the case in which eventually R 3> X s as well. In this situation, 
the increase in coverage is due mainly to the particles deposited within the capture zone of 
the islands, which we assume to be relatively large. 

For times long enough for these conditions to prevail, the evolution of the coverage will 
be given by 

-~(X S + X* S )RFN, (52) 
and the number of islands will evolve as 

^ ~ FVeX| - (53) 
Thus, solving for the coverage, we find that 

6~F A X 2 s (X s + X* s ) 2 T e t\ (54) 

Once again we assume that coalescence becomes important at a value of 6 of the order of 
unity, so that the coalescence time in this regime will be given by 

tc-^XUXs + X^Te)- 1 ' 3 , (55) 

and the maximum number of islands will be given by 

N -'i-0hT- (56) 

Note that ignoring the particles that land on the islands, i.e. taking X s = 0, as if there 
was an infinite edge berrier, yields the same scaling as would be expected if X s 3> X s , and 
also the same scaling as would be expected in homoepitaxy, where Xs = X s . Thus, ignoring 
these " second layer" particles is not crucial for the calculation of the exponents, at least for 
a large range of the values of the parameters. 

If we now assume that eventually R ^> Xs, but that X s ^> R, then once again the 
coverage will evolve as 
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? ~ F9. (57) 

Coalescence then occurs at at a time t c ~ 1/F, and the maximum number of islands at that 
time will now be given by 

N max ~ Fr e X 2 s (58) 
This is the result obtained by Venables, as we discuss later (section VI). 

B. Low evaporation rates 

If r e is further increased, then the system does not reach the regime characterized in the 
previous section. As mentioned above, this regime holds as long as R ^> X$ beyond certain 
point in the evolution of the system. We can estimate the value of r e at which this condition 
never holds, i.e. when X$ ~ R ma x, where the maximum island size R ma x can be estimated 
through 

Rmax ~ (^-) V2 ~ r-Wp-WD-W (59) 



' max , 



Hence, the evaporation time above which the condition R 3> X s never holds is such that 

X, = (Dr e f 2 ~ R max ~ r^p-WD-V* ( 60 ) 

that is, 

rf) ~ (FD 2 r V3 = Zcc 2 r (61) 

For evaporation times above the monomers are expected to travel distances much longer 
than the largest typical island sizes. At this point it is important to realize that our criterion 
for the onset of coalescence, 9 ~ 1, is equivalent to the requirement R ~ £, where £ ~ N^ 1 ^ 2 
is the typical distance between islands. Thus, for evaporation times much larger than , a 
great part of the evolution of the system takes place with the monomers traveling distances 
which are greater than the inter-island distances. During this phase we expect the effects of 
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evaporation to be negligible. In effect, if r e > rj 2 ), the kinetics of the monomer density is 
no longer determined solely by the evaporation, but rather, it is eventually determined by 
the aggregation processes. 

At very short times the number of islands is expected to be very small. Therefore, in the 
early stages of evolution, they cannot affect the monomer density. Thus, the early stages of 
evolution are expected to be similar to those of the previous regime, for example: 

p~Fr e , N~F 2 X 2 s r e t 6~F 3 Xlr e t 2 . (62) 

This situation is expected to hold until the number of islands is large enough, so that the 
term <JiN is no longer negligible when compared to F. That is, until a time t xx such that 

a t N ~ DF 3 X 2 r 2 t xx ~ F, (63) 

from which we obtain 

t xx ~ (F 2 Xlr e )-\ (64) 

As expected, the number of islands at t xx is N{t xx ) ~ l/ATf, indicating that evaporation 
effects "turn" off once the typical inter-island distance is smaller than X$- For times beyond 
t xx all evaporation effects become negligible and our rate equations reduce to those given 
by Tang || for a system without evaporation (plus the direct impingement terms that are 
negligible in this limit) for which the results are well known. Namely the time at which 
coalescence occurs is given by 

t e ~ l/F (65) 

and the maximum number of islands nucleated on the surface is 

N max ~ ( -J (66) 

We can now summarize now the various regimes obtained. There are three principal 
regimes which are spanned as the evaporation time r e decreases. We call them, in this 
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order, the complete condensation regime where evaporation is not important, the diffusion 
regime where islands grow mainly by diffusive capture of monomers and finally the direct 
impingement regime where evaporation is so important that islands can grow only by cap- 
turing monomers directly from the vapor. Within each of these regimes, there are several 
subregimes characterized by the value of X s , the desorption length on top of the islands. 
We use ice = {Ft)~ 1 ^, the typical island-island distance when there is no evaporation and 
Rmax as the maximum island radius, reached at the onset of coalescence. 
complete condensation Xs ice 

N max ~ F^r 1 / 3 for any X* (67) 



diffusive growth 1 <C X s <C ice 

<{FXlT e yl*{X s + X* s )-W if X s <^R max (a) 



N m ax ~ < 



Fr e X 2 s 



if X* s ^>R max (b) 



(68) 



with R max - 



(Xs + X s ) 1 ^(FX s T e ) - 1 / 3 , which gives for the crossover between regimes 
(a) and (b) : X* s (crossover) ~ (FXfr e )~ 1/2 . 



direct impingement growth X s <C 1 

f(Fr e ) 2 / 3 if X^«l (a) 



N„ 



(Fr e ) 2 / 3 X*- 2 / 3 if l«X*«it> ma;E (b) 



(69) 



if X*^>R rnax (c) 



with R 

max rs "' 

(Fr e )- 1 / 3 X^ 1/3 , which gives for the crossover between regimes (a) and (b) 
X* s (crossover) ~ (Fr e ) -1 / 2 . 
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We note that these equations agree with the scaling analysis presented above (Equations 
[| and 0) in the appropriate subregimes (Xg < X s ) ■ 

V. COMPUTER SIMULATIONS 

In the following paragraphs, we test the assumptions and predictions of the analysis given 
in the preceding sections in the special case i* — 1 and no contribution from atoms deposited 
on top of the islands. As we stressed above, the exponents observed without contribution 
from the second layer are the same as those observed for a large range of parameters. We 
also show results that are not attainable from this mean-field calculations, namely the island 
size distributions. 

Our computer simulations generate sub-monolayer structures using the four processes 
included in our model (see the introduction). Here we take r = 1 as the time scale of our 
problem. The monomer diffusion coefficient is then given by D — 1/4. We use triangular 
lattices (six directions for diffusion) of sizes up to 2000 x 2000 with periodic boundary 
conditions to avoid finite size effects. For simplicity, in these simulations, the atoms deposited 
on top of existing islands are not allowed to "fall" on the substrate. We have checked that 
allowing the atoms to fall down the substrate does not change significantly the results. 

The program actually consists of a repeated loop. At each loop, we calculate two quan- 
tities p drop = F/(F + p{\ + i)) and p dif = (p/r)/(F + p(± + ±)) that give the respective 
probabilities of the three different processes which could happen : depositing a particle (de- 
position), moving a particle (diffusion) or removing a particle from the surface (evaporation). 
More precisely, at each loop we throw a random number p (0 < p < 1) and compare it to 
Pdrop and pdif- If p < Pdrop, we deposit a particle; if p > pdrop + Pdif, we remove a monomer, 
otherwise we just move a randomly chosen monomer. After each of these possibilities, we 
check whether an aggregation has taken place and go to the next loop (for more details, see 

§)• 
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A. Checking the assumption : constant coverage for saturation 



A major assumption made in the theoretical treatment is that the maximum of the 
island density is reached when the coverage attains a constant value. We note that this 
assumption is equivalent to the one used in the scaling analysis, since, for compact islands, 
9 ~ NttR 2 ~ (R/l) 2 . It is essential then to check this assumption first. Figs |3]a-c show that 
this it is justified, since 9 max , the coverage at which the maximum island density is reached, 
does not vary systematically with F, r e or e. e is defined as (1 + Xs)Xg (j^j and indicates 
the importance of evaporation : e < 1 means that evaporation is significant, while e ^> 1 
indicates that we are in the regime of complete condensation (see Eq. [2~3|). 
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FIG. 3. Coverage (6 max) a t which the island density reaches its maximum. The different 
figures show the (non) evolution of max versus (a) the flux F (b) the evaporation time r e and (c) 
the evaporation parameter e. The solid lines indicate our assumption, namely that the maximum 
island density is reached at a constant coverage. The dashed lines show linear fits of the data. In 
all these graphs, there is no significative evolution of the coverage at saturation as a function of the 
deposition parameters, as can be seen in the coefficients of the fits : (a) max = 0.13— 0.00085Zog(-F); 
(b) m ax = 0.15 — 0.00069/o#(T e ); (c) 6 max = 0.14 — 0.0012/og(e). We only notice an increase of 
the fluctuations when the evaporation becomes important since the systems contains only a few 
islands. These fluctuations are due to the fact that each point in the figures represents a single 
run, with no averaging. 

In the non-evaporation case, it has also been recognized that 9 max is independent on the 
flux or diffusion rate Actually, Villain and co-workers || use this criterion to find the 

flux and diffusion dependence of N max . We think that the combination of two effects can 
lead to a constant 6 max . First, at coverages close to .15, islands occupy enough surface to 
capture rapidly the landing monomers, which prevents nucleation of new islands. Second, 
at this coverage, islands begin to touch and coalesce (R ~ £), thus starting the decrease in 
island density. We note that 9 max is not constant when large islands are allowed to move, 
even in the non-evaporation case [Q. 

Having confirmed our main assumption, we now turn on to the evolution of the maximum 
island density N max as a function of the different parameters F and r e (remember that we 
take t = 1 as the time scale of our problem). 
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B. Checking the crossover scaling 



Before looking in detail into the different regimes predicted by Equations E% B8| and 



69, we summarize our simulation results in Fig. |j. We show there all our data (more than 
200 points) for N max as a function of the parameters. Our scaling analysis predicts that the 
data should fall into a single curve, given by Equation We see that the data remarkably 
confirms our analysis, over more than 20 orders of magnitude. This gives us confidence on 
our entire approach and its predicted exponents, which we now turn on to check in more 
detail. 
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FIG. 4. Universal function rescaling all our data. As predicted by Equation |23|, the nor- 
malized island density N max Xg follows a single curve as a function of the evaporation parameter 
e = %X S 5 {1 + X S ). The solid curve shows the function predicted in the text (Equation ^4|), while 
the circles represent the results of the computer simulations. 



C. Checking the exponents 



The object of this section is to check that the results summarized in Equations ^7\, 
and B91 are correct. 
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1. Scaling of the maximum island density as a function of incident flux 

Figure |5| shows the evolution of the maximum island density as a function of the flux 
for different evaporation times. Each of these curves is different from the others, since they 
correspond to different evaporation times. However, according to our preceding analysis, 
they should all present a transition from the low evaporation regime to the high evaporation 
regime. This can be detected by a change of slope, from N max ~ F 2 / 3 in the high evaporation 
regime to N max ~ F 1 ' 3 in the low one (Eq. |rJ8| a and |57|) . Of course, this regime change 
does not occur for all the curves at the same value of the Flux, since the parameter that 
determines that change is not the Flux but rather X| = r e / r. Figure [| shows our predictions 
are correct, at least concerning the Flux evolution of the maximum island density. We now 
turn to the other variable, the evaporation time. 




log (Flux) 

FIG. 5. Evolution of the maximum island density as a function of the flux for different 
evaporation times. The number next to each symbol corresponds to the log(r e ) value for that set 
(remember that r = 1). The solid lines show the expected values for the exponents : 2/3 for low 
values of the flux (evaporation is significant) and 1/3 for higher values (complete condensation). 
The dashed line shows the prediction of other authors (N max ~ F x l 2 for low fluxes). See the 
discussion (section VI) for details. 
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2. Maximum island density as a function of evaporation time 



We show in figure ^| the dependence of the maximum island density on r e . We notice that 
for high enough evaporation times, the island density tends to become roughly constant, as 
predicted by our calculations. For lower values of r e , N max changes rapidly. Our analysis 
predicts two regimes : for 1 r e F~ 1//3 , we expect N max ~ r e , (Eq. |B£|a) while for 
r e < 1, we expect N max ~ r e 2//3 (Eq. |rJ§|a). This last regime is clearly seen for the curve 
obtained for a flux F = 10~ 6 (squares, solid line). The intermediate regime is difficult to see 
because of the crossovers with the two other regimes. However, taking a very low value for 
the flux (F = 10" 11 , filled triangles), we can see that the slope in this intermediate regime 
is close to 1 (dashed line). 
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FIG. 6. Maximum island density as a function of the evaporation time for different fluxes. 
The number next to each symbol corresponds to the log(l/F) value for that set. The solid line 
shows the expected value for the exponent : 2/3 for low values of r e (evaporation is significant). 
The dashed line shows the exponent expected for the intermediate regime. 
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3. Direct impingement regime 

We have also checked the exponents obtained when X$ «C 1, in the direct impingement 
regime. Eq. |69| a predicts that the maximum island density scales with the product Fr e 
with an exponent 2/3. This is confirmed by our computer simulations (Figure 0), over more 
than 6 orders of magnitude. 
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FIG. 7. Dependence of the maximum island density N max on Fr e in the direct impingement 
regime {Xg -C 1). The solid line shows our prediction (N max ~ (FY e ) 2 / 3 for extremely low 
evaporation times, Eq. |69|a). We see that our prediction remarquably fits the data over more than 
6 orders of magnitude. The numbers indicate the values of logarithms of the evaporation time r e 
or the flux F for the different sets of data. 




4- Condensation coefficient at the maximum island density 

A last test for the analysis is presented in Fig ^. We show that the dependence of the 
sticking coefficient S at the saturation island density (5* = 9 max /Ft c ) follows Eqs. ER 03 



and |65]. By constructing an evaporation parameter rj = FX|(1 +X|), we can group all the 
regimes in a single curve : for the complete condensation regime, S ~ i]° ~ 1, for the others 

S ~ 7] 1 / 3 . 
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FIG. 8. Evolution of the sticking coefficient at the maximum island density as a function of 
the evaporation parameter FXg(l + Xg). The solid lines represent our predictions (Eqs. |55| and 
65). See the text for details. 



D. Island size distributions 

Island size distributions have proven very useful as a tool for experimentalists to distin- 
guish between different growth mechanisms |25| , |32"ll . By size of an island, we mean its total 
number of monomers or mass. Unfortunately this information is beyond the reach of the 
simple mean field rate equation analysis presented above || . Nevertheless, the distributions 
can be obtained from the simulations. Fig |] shows the evolution of the rescaled ||7|,|32| island 



size distributions as a function of the evaporation parameter e = (1 + Xs)X s . It is clear 
that the distributions are significantly affected by the evaporation, smaller islands becoming 
more numerous when evaporation increases. This trend can be qualitatively understood by 
noting that new islands are created continuously when evaporation is present, while nucle- 
ation rapidly becomes negligible in the complete condensation regime. The reason is that 
islands are created (spatially) homogeneously in the last case, because the positions of the 
islands are correlated (through monomer diffusion), leaving virtually no room for further 
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nucleation once a small portion of the surface is covered (6 ~ 0.05). In the limit of strong 
evaporation, islands are nucleated randomly on the surface, the fluctuations leaving large 
regions of the surface uncovered. These large regions can host new islands even for relatively 
large coverages, which explains that there is a large proportion of small (s < s m ) islands in 
this regime. 
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FIG. 9. Rescaled island size distributions for different values of the evaporation parameter e 
and different coverages 9 between .05 and .2. p(s/s m ) represents the probability to find an island 
containing s monomers when the average size is s m — /N . If ti s is the number of islands containiug 
s monomers, we have p(s/s m ) = s m *n s /Nt where Nt is the total number of islands. The solid line 
shows the size distribution obtained without evaporation. The first number next to each symbol 
corresponds to log e. The two other figures correspond to log F and log r e respectively (r = 1). 

VI. DISCUSSION 



Other authors have analyzed similar mean-field rate equations to find the growth dy- 
namics and maximum island density in the presence of evaporation [|T5| , p!7H l9|1 . In what 
follows, we discuss the relationship between our work and the preceding studies. In fact, a 
distinction must be made between the result of Stoyanov and Kaschiev, and that of Venables 
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et al. 

Stoyanov and Kaschiev consider that atoms may diffuse on top of an island, with the same 
diffusion length as on the substrate (Xg = Xg in our notations). This corresponds to the 
regimes described by Equations ^ and |68| a (they do not consider the direct impingement 
regime). However, their result in the regime where evaporation is important is N max ~ 
(FTe) 1 / 2 , different from our predictions (see both Eq. |6^ a and Eq. |69|b), which are clearly 
supported by our simulations (see for example Fig. ||). As shown in the following, we think 
that the difference with our result stems from overlooking the fact that only the adatoms 
falling at a distance X s from the island edge contribute to the rate of growth of the island : 
in other terms, their capture cross section is not appropriate. 

Stoyanov and Kaschiev assume, as we do, that the diffusing adatoms have two possible 
fates: either to evaporate or to be captured by an island. Their lifetime Tj is thus 

— = — + DN , (70) 

and the adatom density reads 

p « Fr { « — — . (71) 

H l + Dr e N K 1 

This is Eq. (6.2) of @. 

At this point, the incorrect assumption is made that the rate of growth of the area 
of an island of radius R reads dR 2 /dt ~ Dp ps FDr e /{l + Dr e N) throughout the strong 
evaporation regime (cf. Sec IV). Writing then 

£> p (**+i) ~ ^iv 2 w DpN 2 , (72) 



one finds p % * ~ N 2 , or 



N 2 (l + X 2 Nf « [^X 2 ) . (73) 



This is Eq. (6.12) of |T9|| , and it leads to N « (F/ D) % * l 2 X % * in the "incomplete condensa- 
tion", X 2 N < 1. 
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On the other hand, Venables et al. actually study (though it may not be easily guessed by 
reading their papers : we would like here to acknowledge an anonymous referee for precious 
help) the regime where all the particles deposited on top of the islands immediately fall on its 
border, thus contributing to increase its radius. This corresponds to the subregime described 



Table 1 of Ref. ||18|| ). This regime applies only if diffusion and/or desorption are different on 
the substrate and on top of islands, which cannot be the case in homoepitaxial situations. 
In heteroepitaxial growth, on the other hand, the substrate and the islands of the first 
monolayer are chemically different, and Venables et al.'s assumption may apply; however, it 
is a special situation, not the general rule. The origin of their "initially incomplete" regime 
is more mysterious. We try in the following to understand its origin. 

We now show that the assumption of Venables et al.'s concerning the atoms falling on 
top of an island, leads to their equation for the maximum island density, Eq. (2.17) of Ref. 
fl"8|| . In our notations the latter reads 



where 9 = R 2 N is the deposited dose (coverage) and 9 its value at coalescence. 

Assume that an island grows by capturing all atoms falling on top of it, plus those 
diffusing to it. Then, the rate of growth of the island area is 



by Eq. |5"E|b. Our prediction agrees with their result in the "Extreme Incomplete" regime (see 



N(l + X 2 s Nf% + X 2 S N) « (F/DfX i 



(74) 



dR 2 



FR 2 + Dp^F9/N + Dp, 



(75) 



dt 



where 9 = R 2 N is the deposited dose, or coverage. 



From Eq. (|71|) and the nucleation rate, one finds 




(76) 



Letting fl7|) equal to (dR 2 /dt)N 2 (cf. Eq. |72|), one has 




(77) 
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where #0 < 1 is the coverage at coalescence. Rewriting yields 

« H 1 + X * N ) + X * N ] N^(9o + X 2 S N)N , (78) 

which gives immediately Eq. flT4"|). 

We now turn on the relation of our crossover scaling (section III-C) to recent work by 



Bales ||34j| . Recently, the crossover scaling of the island density between different growth 



regimes on vicinal surfaces (without evaporation) was addressed by Bales [34J], and by 
Pimpinelli and Peyla In particular, Bales claims that "when atoms desorb from the sur- 
face a crossover scaling form identical to" his results "is satisfied with" the distance between 
surface steps "replaced by the average distance X s a monomer will travel before desorbing" 
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We argue that Bales' claim is not true, as it stands. To see this, we need consider what 
happens on a vicinal substrate. We do this using the argument of Sect. III. Of course, 
competition between two length scales occurs in this case, too, the two lengths being the 
adatom diffusion length before capture by an island, £ s , and the step-step distance I. When 
i s <C £, the diffusing adatoms behave as though they were on a flat substrate, and the island 
density is still given by equation (4). In the opposite case, £ s S> £, the adatom density 
is fixed by capture at steps, so that p can be found by replacing X s by £ in Eq. (5), in 
agreement with Bales' statement: 

p » Ff/D . (79) 

However, since evaporation is now neglected, an island grows by capturing all adatoms 
diffusing to its border. This means the island average radius grows as 



R w NFt , (80) 

and since R ~ £ s eAt — t c , one has 

1 « f| « FfN . (81) 
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In the regime of dominant adatom capture by steps the island density at saturation is 
found from 

Dff+ l « — (82) 

tc 

and thus 

N » (-£ 2 ) . (83) 
A crossover between equation (4) and (|83D is thus expected: it has been actually found 



by Pimpinelli and Peyla ||33|| . 

Bales |33] considered the crossover scaling of the island density at low coverage. His 
result for the regime of dominant adatom capture by steps, Eq. (9) of his paper, is found 
(within logs) by letting t c = t = 9/F and i* — 1 in ( [82] ) above, 

N « ^£ 4 9 . (84) 

Note that this cannot be obtained by replacing X s by I in our result, Eq. (7) above. It 
is however obtained by replacing X s by I in Venables' result. This is because Venables does 
not take into account the crucial phenomenon which occurs during evaporation-namely, that 
only the adatoms falling near or at the step edge can contribute to the island growth. To 
neglect this is correct within Venables' assumptions, as well as on a vicinal surface without 
evaporation, where the latter mechanism is replaced by adatom capture by steps. 



Our crossover scaling does not therefore coincide with Bales', and Bales' results |]3"4" 
however modified do not describe homoepitaxial growth with desorption (Xs = Xg). 



VII. SUMMARY AND PERSPECTIVES 

By combining different mean-field analysis and extensive computer simulations, we have 
shown that the presence of evaporation has important effects on the growth of submonolayer 
films. We have investigated the different regimes that arise when the growth parameters 
are varied, and predicted the behaviour of several experimentally accessible quantities such 
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as the island size distributions, the maximum island density and the time at which this 
maximum is reached. In some cases, by measuring these last two quantities we can infer 
the values of both the evaporation and the diffusion times, which are difficult to obtain 
otherwise. For example, if the experiments are carried in the intermediate regime (and this 
can be checked from the sticking coefficient), by using Equations |58| a and [55], we obtain 
r = F 2 (t c N max f and r e = t c N max 2 . 

We have also shown that our model is more general than previous studies of growth in 
presence of evaporation. 

Future directions of study should include more realistic hypotheses for a direct compari- 
son with experiments [16] . In particular, we are presently extending our analysis to the case 
of three-dimensional growth (i.e. leading to cap-shaped islands) and the influence of defects 
on the surface which could act as nucleation centers. 

We wish to thank M. Meunier and C. Henry for helpful discussions, Ph. Nozieres for a 
critical reading of the manuscript and an anonymous referee for precious comments on the 
connection of our analysis with previous works. 
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